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Abstract 

There exists a simple relationship between a quantum-mechanical bound- 
state wave function and that of nearby scattering states, when the scattering 
energy is extrapolated to that of the bound state. This relationship is demon- 
strated numerically for the case of a spherical well potential and analytically 
for this and other soluble potentials. Provided that the potential is of finite 
range and that the binding is weak, the theorem gives a useful approximation 
for the short-distance behaviour of the scattering wave functions. The connec- 
tion between bound and scattering-state perturbation theory is established in 
this limit. 

I. Introduction 

It has recently been shown there is a simple relationship between the normalisa- 
tions of the scattering and bound-state wave functions when the scattering energy 
is continued to that of the bound state 1 . That the two wave functions become mu- 
tually proportional in this limit is known from the standard text books 2 ' 3 , but it 
was unexpected to find that the relative normalisation depended purely upon the 
binding energy. Though the theorem is valid for arbitrary angular momentum, we 
wish here to illustrate some of the results for S-w&ve scattering. Provided the range 
of the potential is short and the binding is light, it will be seen that the extrapolation 
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theorem actually gives useful approximations for scattering wave functions at low 
energies and short distances. 

The S"-wave Schrodinger equation for the interaction of a particle of mass m in 
a potential V(r) may be written as 

u"(r) - a 2 u a (r) = U(r)u a (r) 
v"(k,r) + k 2 v(k,r) = U(r)v(k,r), (1) 

where U(r) = 2mV(r)/h 2 and the binding and scattering energies are —h 2 a 2 /2m and 
h 2 k 2 /2m respectively The corresponding radial wave functions, u a (r) and v(k,r), 
vanish at r = but satisfy contrasting boundary conditions at large distances. For a 
potential of finite range the bound-state wave function behaves asymptotically like 

u a (r)^N as e~ ar (2) 

and is normalised by the integral condition 

/ [u a (r)} 2 dr = 1 . (3) 
Jo 

On the other hand the scattering solution normalisation is determined by its asymp- 
totic behaviour and, if we take real boundary conditions (standing waves), then 

v(k, r) — > ^ sin(fcr + 5(h)) = \ (e im e tkr - e^V^) (4) 

as r — > oo, where 5(h) is the 5-wave phase shift at wave number k. 

In order that v(k,r) might become an acceptable localisable wave function as 
k — > ia, for which the normalisation condition of Eq.(||) would hold, the S-matrix 
must have a pole at the position of the bound state and can therefore be written in 
this region as 

s W = **» = jvw (^|) ^ = N{k) , (5) 

where N(k) is well behaved in the vicinity of k = ia. 

At the bound-state pole the second term in Eq.@ is eliminated to leave the 
desired asymptotic behaviour 



Though this does not establish the relative normalisation of this and the bound- 
state wave function u a (r), the presence of the square root in the denominator of 
Eq.(|6]) is indicative and this factor is contained within the extrapolation theorem 1 , 
which for potentials of finite range yields 

2a(a 2 + k 2 )Y /2 v(k,r)\ = -u a (r) . (7) 



lim 

k— >ia 

It should be noted that this relationship is independent of the shape of the potential 
and is determined purely by the binding energy. It is now easy to see from Eqs.@ 
and (|7|) that, in terms of the asymptotic normalisation constant N as of Eq.(0), the 
residue at the pole of the S'-matrix in Eq.(j5|) equals —iN 2 s . 

As will be seen from the numerical examples given for the spherical well potential 
discussed in §2, the above extrapolation theorem actually provides useful approx- 
imations for scattering wave functions in terms of that of a nearby bound state. 
Provided that the binding is weak, and that neither k nor r are too large, 

v(k, r) « - [2a{a 2 + k 2 )] ^ u a (r) . (8) 

The ratio of the scattering to the bound-state wave functions is studied analytically 
for several exactly soluble potentials in §3, where more precise statements are made 
on the range of validity of the approximations in Eq.([8|) and general trends in the 
deviations noted. One feature seen for the spherical well potential is that there seems 
to be a region of r where all the wave functions, scaled by the square root factor 
as in Eq.(|8|), come together. A simple understanding of this cross-over effect is also 
attempted in §3. 

The low energy proton-proton system is not quite bound but there is a virtual 
state with small but negative a situated close to physical energies. Though in such 
cases there is no bound-state wave function to set the scale, the energy dependence 
of the scattering wave functions is well approximated by the square-root factor as in 
Eq.(||). This behaviour is clearly demonstrated by the spherical well and it is seen 
from these examples that it is the nearest singularity which is dominant. 

Though the long range of the Coulomb potential violates the assumptions made 
when proving the extrapolation theorem of Eq.(M), it nevertheless remains valid for 



all the Coulomb bound states, as shown in §4. However the accumulation of bound 
states in the vicinity of zero energy means that Eq.@ is of no use in representing 
the scattering functions even at low energies. There is no nearby dominant pole. 



processes at large momentum transfers, where the behaviour of the wave function 
at small distances is important 4 ' 5,6,7 . For example it allows one to estimate the cross 
section for pp — > pmr + in terms of that for pp — > dir + , where the final neutron-proton 
triplet pair are fused to become a deuteron. Another simple example discussed in §5 
involves the relationship between perturbation theory applied to bound and scatter- 
ing state problems. A consistent picture is found when using Eq.(|||) in the presence 
of a lightly bound state. A summary of our main conclusions is to be found in §6. 

II. Numerical investigation for the spherical well potential 

As a first illustration of our ideas, consider the spherical well potential 



For r > a the scattering and bound-state wave functions are proportional to 
sin(fcr + 5) and e~ ar respectively, whereas inside the well they are rather sin(/tr), 
where k 2 = Uq — a 2 or h 2 = Uq + k 2 . It is trivial to match boundary conditions at 
r = a to obtain an implicit equation for the energy eigenvalues and hence the wave 
functions 2 . 

The approximation of Eq.(^) is expected to work best for a loosely bound state, 
and we start by treating the case where there is just a single Is state. Taking units 
where the potential radius a = 1, a potential strength of Uq = 2.8 yields a = 0.159 
and the normalised full bound-state wave function ip a {r) = u a (r)/r is plotted in 
fig. 1.1. Also shown there are the corresponding scattering wave functions, modified 
by the square-root factor of Eq. (|8|) , 



Eq.(^) is of great help in estimating the effects of final state interactions in nuclear 




r < a , 
r > a . 



(9) 




1/2 



v{k, r)/r , 



(10) 



at wave numbers of k = 0.1, 0.2, 0.5, and 1.0. 
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At r = the agreement between the various functions is very good for low values 
of k, though the scattering ones lie above that of the bound state and steadily increase 
with k. Though the energy dependence of the functions in this region seems to be 
roughly k 1 ^ 2 , too much should not be read into this. The scattering functions must of 
course be orthogonal to the bound state and this implies that they start to oscillate 
before the bound state has died out. This is seen clearly in the figure, as is the fact 
that the oscillations set in earlier for higher fc-values. What is less evident is why 
all the curves seem to come close together for r « 0.66. This has the practical effect 
however of limiting the deviations between the scattering and bound-state functions 
over the range of the potential. 

If the potential strength is increased in order to generate a deeply bound Is state 
in addition to a 2s state with the same binding as in fig. 1.1, then the latter has the 
node demonstrated in fig. 1.2. The scattering functions defined by Eq.(^[) still fall 
very close to that of the bound state with the crucial difference that they now lie 
below for small r. A cross-over region still exists but at higher values of r. 

Fig. 1.3 shows that if the attraction is reduced a little such that the 2s state 
becomes a virtual one with negative a (a = —0.159) then, as expected, the scattering 
wave functions still retain their dominantly 2s character, though there is then no 
bound state to set the scale and r-dependence. The wave function of such a virtual 
state increases like e + ' a ' r at large distances and hence is non-normalisable. It is 
therefore clear that it is the nearest singularity of the scattering amplitude which 
governs the behaviour of the wave function at low k. Though Eq.(|7]) is valid for 
all the bound states of a problem when extrapolating to the pole, Eq.(^) is only a 
useful representation for a nearby singularity. An important practical case of this 
kind is the proton-proton system where the (just) unbound virtual state dominates 
the low energy wave function and allows us to derive the energy dependence of the 
pp — > ppn° cross section near threshold 5 . 
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III. Analytic expressions for soluble potentials 



In addition to the spherical well potential, there are several other potentials for 
which exact solutions to the Schrodinger equation can be obtained. In such cases 
it is interesting to study the deviations from the extrapolation theorem analytically 
in going away from the bound-state pole at k — ia to see the dependence upon the 
parameters. Define therefore the ratio 

r i 1/2 

R(k,r) = - [2a{a 2 + k 2 )\ v(k,r)/u a (r) (11) 

In order to satisfy the extrapolation theorem R(ia, r) = 1 and, provided that the 
potential is of finite range, the ratio function can be expanded as a power series in 
a 2 + k 2 . 

oo 

R(k, r) = 1 + £ Mr) (« 2 + k 2 ) n . (12) 

n=l 

In general this will have a finite radius of convergence since R(k, r) diverges at the 
position of any other bound state. 

A. The Yamaguchi potential 

The strength A of the non-local but separable Yamaguchi potential, 

nry) = -*(^P) (^) , as) 

can be adjusted to give a solitary bound state at E — —h 2 a 2 /2m. 
A straightforward calculation then shows that to second order in r 



B. The Bargmann potential 

The Bargmann potential defined by 

U(r) = V(r) = -2(P 2 - a 2 ) ^cosh(/?r) - ^ sinh(/3r)^ , (15) 

with (3 > a > 0, has precisely one bound state with E = — % 2 a 2 /2m. 
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Taking the ratio of the scattering to bound-state wave functions at low energies 
and distances then gives 



C. The spherical well potential 

The solutions for the spherical well potential U(r) = —Uo9(a — r), discussed 
numerically in §2, also lead to analytic expressions for the ratio function. Thus 

„ . . a 2 (4 + ad) 1 , / A \ , 

Rl r = 77TT — 7 - 7777 k ~ r r + °( r ' 17 

4(1 + aa) 4(t7o — ol ) o v ' 

where the binding energy E = —h 2 a 2 /2m is determined from the implicit equation 

(Uq - a 2 ) 112 cot (a(U - a 2 ) 1/2 ) = -a . (18) 



D. Conclusions 

It is clear from the examples given above that for a lightly bound state (a 
£7(0)), the expansion parameter at the origin is proportional to (a 2 + k 2 ) times 
the square of the range of the potential, though the coefficient depends upon the 
shape of the potential. However the r-dependence of the expansion coefficients R n {f) 
can be obtained directly by expanding the potential in powers of r and integrating 
the Schrodinger equation starting from r = 0. It is therefore no accident that the 
coefficient of r 2 in Ri(r) is the same for the Bargmann case of Eq. ([!(]) as for the 
square well of Eq. (|TT|) . This result is independent of the potential provided that this 
is finite at the origin. In the Yamaguchi example of Eq.(jl4|) the odd powers of r 
arise from the divergence of the potential at short distances. 

It was noted in §2 that in the case of the first lightly bound state in a spherical 
well potential the scattering function lay above the bound-state function, whereas if 
it were the second bound state then the scattering function lay below. This change 
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in sign of i?i(0) in Eq.flTTD arises because in the first case Uo ~ 7r 2 /4a 2 , whereas in 
the second U m 9n 2 /4a 2 . 

In the numerical examples shown in §2, where the deviations from the extrapola- 
tion theorem are not too large, a cross-over region was noted where all the functions 
are very close. That the scattering function cross the bound state, in a region where 
the latter wave function is large, follows immediately from the condition that the two 
be orthogonal when integrated over r. This position will be stable with respect to 
changes in k provided that the first term in the power series expansion in (a 2 + k 2 ) in 
Eq. (|i~2"D is dominant and in such an event the cross-over point would be determined 
by the condition that Ri(r) = 0. For (3 = 1.0 and a = 0.1 in the Bargmann case 
one would expect that at low energies the cross-over would occur when r « 1.74 
if we keep only the quadratic term in r 2 in Eq . (|I~6"1) , but r « 1.68 with the quartic 
term. These are to be compared to the numerical value of r ~ 1.522. Of course the 
cross-over point should move to smaller values of r as the energy is increased since 
the scattering wave functions start to oscillate faster and the node in r is shifted to 
the left. 

The above arguments do not provide any explanation as to the sign of R\(0) and 
for this the spherical well case is particularly illuminating. In the case of fig. 1, where 
one has only one (lightly)-bound state, then i?i(0) is positive. This is also true for 
the other two potentials discussed in this section. However when the depth of the 
well is increased so as to make the lightly bound state the second one then -Ri(O) 
changes sign. This is in the right direction to make the extrapolation theorem of 
Eq.(|7|) be valid also for the deeply bound state at more negative values of k 2 where 
R(k, r) has to diverge to +oo. On the other hand, when there is only one bound 
state the only singularities that R(k, r) can be simulating are those associated with 
the potential itself. 
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IV. The Coulomb potential 



Due to the long range of the Coulomb potential the infinite number of bound- 
state eigenvalues have an accumulation point at energy zero. Nevertheless it is easy 
to see explicitly that the extrapolation theorem is still valid for all the bound states. 

Define full scattering wave functions as ijj(k,r) = v(k,r)/r, and similarly for the 
bound state. The value of the square of the real Coulomb wave function at the origin 
in an 5-wave scattering state is just the Gamow factor 

^ °)] 2 = ^TT7 = ^ t- 1 + COth(7T77)] , (19) 
where in the attractive case the Coulomb parameter r] = —me 2 /h 2 k. 



The infinite number of bound states correspond to the poles of Eq.(19) when 
77 = in, where n a non- vanishing integer 8 . These can be made explicit by recasting 
the equation as 



ty(k, 0)] 2 = -7T77 + 1 + 2 J2 ; 1 „ 2 —1 ■ (20) 

71=1 L 



1 

+ n 2 /f] 2 

If we write the virtual momenta of the bound states as a n = me 2 /nh 2 , it is well 
known that the value of the n'th bound-state wave function at the origin is 



[Va0)] 2 = 4c4, (21) 



so that Eq.(EDI) may be written as 



Our theorem of Eq.fl7|) is still formally valid since, in the vicinity of any of the 
Coulomb bound states, just one term in the sum will be important and for this the 
residue is clearly correct. Despite this, it is of no practical use in representing the 
scattering wave function at low energies since no single pole will then dominate and, 
as can be seen from Eq . fl22|) , the Coulomb potential produces an extra non-trivial 
term singular at k = 0. Thus the expansion of R(k,r) in powers of (a 2 + k 2 ) as in 
Eq.(|T2p is not valid. 
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V. The two-potential formula 



We here show the relation of our approximation to some of the results of standard 
perturbation theory. Denote by /o the S-wave scattering amplitude corresponding 
to a potential Vo- If a small extra potential V\ is added, the change in the amplitude 
to first order in V± is given by the two-potential formula 3 

A/«-^e«*(x*l^i|x*>- (23) 
n 

The factor involving the phase shift S for the potential V is introduced because Xk 
is the real scattering solution. Note that the integration in the expectation value is 
only over the radial coordinate r. 

Assume now the potential V\ to be of short range such that the conditions leading 
to Eq.flj]) are applicable. If the potential V has one weakly bound state Xb, then Xk 
is related to it through the approximate relation of Eq. (H) and this leads to 

A/ W - e ™° i ^\VAXb) = ~e 2 ^ ■ (24) 

In the above we have used perturbation theory to identify AE as the change in the 
binding energy E = —a 2 l /2mh 2 of the bound state to first order in V\. In this limit 

2ikmAE 
a (a.Q + k 2 )h 2 

For a loosely bound state, the S-matrix is given by the scattering length approx- 
imation 

e^ = '^±% and e 2 ^ = ^ • (26) 
ao + ik a + ik 

The correction in Eq . (|25|) is valid only to leading order in AE and introducing 
this into Eq.(p6|) and keeping only such terms, we find 

mAE s 



e 2iS w 1 + 2ik(f + Af) = e 2iSo 



1 



(25) 



a = a [ 1 - (27) 



and a binding energy of 

E = — — — = E + AE , 28 
2m 

as expected. This explicitly shows the consistency of first order perturbation theory 
for bound-state energies and scattering matrices with our approximate wave function 
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in the case of a weakly bound state. 

VI. Summary and conclusions 

We have shown for several soluble examples that it is possible to approximate 
the real scattering wave function quantitatively in terms of that of a nearby bound 
state. This goes further than the theorem of Eq.(^) in that it shows that it is a 
robust extrapolation provided that the binding is weak and that the wave number 
k and distance r are both small. Though we have worked entirely with real wave 
functions, this is sufficient for applications in final-state-interaction theory . 

Provided that the bound-state pole is close, one can use Eq.@ to obtain the com- 
plex scattering wave function with outgoing boundary conditions in the scattering 
length approximation 

^+\k, r) = ^j{k, r) e iS{k) « -(2a)~ 1/2 ^ a {r)/{a + life) . (29) 

Though we have a qualitative explanation for the deviations from the approxi- 
mate relationship of Eq.@, its magnitude at r = must depend upon the structure 
of the potential and the singularities of the scattering amplitude. Fortunately the 
cross-over phenomenon, whereby all the functions come very close at some region 
within the range of the potential, limits the deviations significantly. As explained in 
§3, this cross-over phenomenon is very general and exists also for realistic nucleon- 
nucleon potentials, such as that of the Paris group 4,9 . It is also demonstrated by 
numerical resolution of the Schrodinger equation for various finite-range potentials. 
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Figure 1: Bound-state (solid line) and scattering wave functions, modified as in 
Eq. fllCf) , for a spherical well of radius 1.0. The scattering states are evaluated at 
k = 0.1 (dot-dashed line), k = 0.2 (closely spaced dots), k = 0.5 (dashes), and 
k = 1.0 (widely spaced dots). (1) For a depth of Uo = 2.8, the resulting Is bound 
state has a = 0.159, and the deviations at r = are positive and increase with 
k. (2) A potential of strength Uq = 22.547 has a 2s level with the same value of 
a = 0.159. The scattering functions now lie below that of the bound state at r = 0. 
(3) Reducing the potential strength to Uq = 21.913 turns the 2s level into a virtual 
state with a = —0.159. Nevertheless the scattering functions retain their dominantly 
2s character at short distances and illustrate that it is the nearby singularity, bound 
or virtual, which governs the energy dependence of the scattering wave functions. 
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